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ABSTRACT 

Using the Star Track binary population synthesis code (jBelczvnski et alj |2008[ ) we 
model the population of binary pulsars in the Galaxy. We include the detailed treat- 
ment of the spin evolution of each pulsar, due to such processes as spin-down, spin- up 
in the accretion events and magnetic field decay. We also model the spatial distribution 
of the binary pulsars due to their propagation in the Galactic gravitational potential. 
Newly born neutron stars are given natal kicks. This synthetic pulsar population is 
compared to the observed sample of binary pulsars taking into account the selection 
effects of detection in the radio band, to determine the best fit evolutionary param- 
eters. With these parameters we determine the properties of the double neutron star 
binaries detectable in gravitational waves by the high frequency interferometer LIGO 
and VIRGO. In particular, we discuss the distributions of chirp masses and mass ratios 
in the radio and in the gravitational wave selected sample. 

Key words: binaries: general - stars: neutron - stars: statistics. 



1 INTRODUCTION 

Double neutron stars (DNS) are remarkable objects in as- 
trophysics and their studies uniquely contribute to a number 
of areas. The monitoring of PSR 1913-1-16 provided strong 
limit s on the General Theory of Relativity (|Hulse fc Tavloij 
1 19751 ). The recent discovery of the first system with two 
observable pulsars J0737-3039 (|Burgav et al.l (|2003l ) and 
iLvne et al. has been a breakthrough in pulsar stud- 

ies. The eclipses in this system provide insight in the p hysics 
of the magnetosphere ((Lvutikov fc Thompson! 120051 ) ■ The 
possibility of measurement of the moment of inertia of pul- 
sar A in this system may lead to very strong constraints on 
the neutron star equation of state. DNSs are the only com- 
pact object binaries that we can study and deduce from their 
properties that there is a population of merging binaries - 
sources of g ravitational waves f or the LIG O and VIRGO ob- 
servatories lAb ramovici et al.1 1 19921 ) and iBradaschia et all 
l[l993)) 

The origin of DNSs has be en first worked out by 
iFlannerv fc van den Heuvell l|l975h in the context of the the 
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Hulse Taylor pulsar PSR 1913-1-16. T he birth rate and prop- 
erties of DNSs wer e inves ti gated by iPiranI (Il992l). a nd by 
lArzoumanian et all (Il999l ). iBelczvhski fc Kalogeral (120011 ) 
pointed out that there may exist a large population of ultra- 
compact DNS binaries formed with an additional common 
envelope stage, and not detectable in the radio. The exis- 
tence of this population relies strongly on the result of the 
common envelope phase with a star on the Hertzsprung gap. 



The inve stigation of pulsa 
tempted by iGunn fc Ostriken 



evolu tion was first at- 

1970h soon after the 

papers concentrated on the 

sample of pulsars with the 



discovery of pulsars. Most 
modelling of the observed 
inclusion of all selection effects and modeling of the pulsar 
evolutionary parameters such as initial distri bution s , ma] 

(■Large (l97 
(1985) , Stollma 



netic field d ecay, evolution of inclinatio n etc 
iTavlor fc Manchester ( 1977 



. | Lyne et al. , _ _ _ , _ 
1987), Emme ring fc CheTOheTll 19891). iBhattacharva et al.l 



(il992). iHartman et al 



Arzoumanian et al. 



Gonthier et al] (|2004l ) 



1 20021) ■ 



Kalogera et al.l 



Gonthier et al.l 



(120011 ) 
ll2002h 

Faucher-Giguere fc Kaspil (|2006| ) 
). These efforts were mostly concentrated on the classical 
pulsars and recently IStorv et al.l |20o3) have presented 
models of the population of recycled millisecond pulsars. 
However, they do not take into account the effects of binary 
evolution leading to formation of these objects but assume 
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initial injected distributions of magnetic fields and spi n 
periods. For a review on binary pulsars see iLorimerl l|2008h . 

The evolution of binaries leading to formation of double 
compact objects has b ee n inv est igated by many authors 
jGnusareva fc Lipunov 1 19851: iBrandt &: Podsiadlowskil 



10' 



19951: IJorgensen et al] Il995l: 



199S 



Lipunov et al 



Popov et al.l 



200C 



Portegies Zwart fc YungelsonI ._ 

Dewi et al ll2005l : iBogomazov et al.ll2005l : iDewi et al.ll200e 
Kiel et al.ll200i ) 



1995 



Here we will be using the Star Track binary pop- 
ulation synthesis code which was presented and used in 
a series of papers by Belczynski and collabo rators, see 
iBelczvnski et all (|2002l ). iBelczvnski et all l|2008l ') and refer- 
ences therein. The code has been used to investigate the 
possible binary progenitors of gamma ray bursts, and also to 
the rates and properties of binary gravitational wave sou rces 



BuUk fc Belczvriski 2003; Bulik ct a l. 2004). Bulik ct al 



20041 ) presented a calculation of the expected mass spec- 



trum of merging compact object binaries including the dou- 
ble neutron stars. They hinted at a possibility of detecting a 
large number of non equal mass NS-NS binaries, and showed 
that they may not show up in the ra dio sample due to se- 
lection effects l|Belczvnski et aLll2002l ) . 

In this paper we present the merger of two models: 
binary population synthesis and pulsar evolution models 
and compare them with the available data. In section 2 we 
present the model of the binary evolution and the pulsar evo- 
lution we use. Section 3 contains the results and comparison 
with the observations. Section 4 is devoted to comparison 
of the radio and gravitational wave selection effects, while 
Section 5 contains the conclusions. 



1.1 Properties of known binary pulsars 

There are currently 10 known pulsars in nine NS-NS bi- 
naries. In one case, both stars in the binary are visible as 
radio pulsars (J0737-3039). At Fig. [1] we present the spin 
period - spin period derivative diagram {P — P) for this bi- 
naries. Most of the objects are concentrated in the millisec- 
ond pulsar region. However there are two noticeable outliers: 
J07373-3039B, the companion of the J07373-3039A in the bi- 
nary pulsars system; and J 1906-0746, a young pulsar with a 
likely neutron star companion. Detailed prope rties of those 
system s ca n be found in T ab. [T]([st airs (2004) F aulkner et al.l 
l|2005l ) and iKasianI (|2008l )). The binaries in Tab. [1] are or- 
dered according to the merger time (time remaining to the 
coalescence). Three binaries at the bottom of the table have 
merger times much longer than 10 Gyr. We list the spin pe- 
riod and its derivative, the masses of the neutron star, as 
well as the present orbital parameters. The distribution of 
the orbital periods spreads over two orders of magnitude, 
with no clear evidence of clustering in these range, and the 
orbits of all systems are significantly eccentric. The masses 
of the seven objects with the merger times below 10 Gyrs 
are determined very well and are in the range between 1.25 
and 1.44 Mq. The masses of the neutron star in the remain- 
ing three binaries are not so well constrained and may even 
lie outside this range. 
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Figure 1. Observed binary pulsars on P — P diagram. Name of 
each pulsar is given next to it. We currently know about 10 pulsars 
in NS-NS binaries. In one case both are visible as radio pulsars 
(J0737-3039). The death lines are described in section 2.3.8. 



1.2 Radio pulsars and population observable in 
gravitational waves 

Due to the selection effects observed sample of NS-NS bina- 
ries might differ from the intrinsic population. There is also 
a possibility that the population observed in radio will have 
different properties that the one observed in gravitational 
waves. The most obvious selection effect in radio is drop of 
observed flux proportionally to distance squared from the 
observer to the pulsar. Our radio-telescopes have limited 
sensitivity so we will not see the weakest and/or farthest 
pulsars. The binary pulsar population observed in radio is 
restricted so far to our Galaxy. 

In this paper we assume that all pulsars with observed 
flux below ~ 1 mjy at 400 MHz are not detectable in ra- 
dio. We also try to incorporate in our calculations a second 
important selection effect directly connected with the fact 
that the pulsar s considered are in binaries. According to 
iFaulknej l|2004l ). while searching pulsars with stack search 
and phase modulation methods it is difficult to detect radio 
emission from neutron stars which are in binaries with or- 
bital period within range [0.3 h;4 h]. In general it is due to 
the fact that signal-to-noise ratio drops drastically due to 
variation of pulsar period in during the observation. There 
are of course different methods of detecting radio pulsars 
which are sensitive for pulsars of such properties. Further- 
more this orbital period range was calculated for a pulsar 
with the spin period P = 9 ms l|FaulkneJl2003 ). and it may 
vary as a function of P. 

However, our aim is to see how this selection effect influ- 
ences the observed properties of the population. Binaries de- 
tectable in gravitational waves will be those with relatively 
small orbital separation. Namely, we assume that binaries 
with coalescence time longer than 10^° years a r e not visi- 
ble in gravitational waves. In both iBulik et al.l l|2004h and 
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Table 1. Properties of the observed pulsars in NS-NS binaries. The table contains names, spin period, spin period derivative, 
orbital period, mass of observed neutron star, mass of the companion, eccentricity of the orbit and time to merger. All given digits 
are significant. Errors, where given, are I — a errors. 



Name P [ms] P [ss-i/lO-l**] Port M M^t, [Mq] M^^p [Mq] e Tmrg [Gyr] 

J0737-3039A 22.70 1.74 2.454 l-337j;o:oo5 l-250tooog 0-088 0-085 

J0737-3039B 2773 8.8 x 10^ 2.454 1-2501q l-337+°;Qog 0.088 0.085 

B2127+11C 30.53 4.99 8.05 l-35lg;°* l-36lo.o4 0.681 0.2 

J1906+0746 144.07 2.028 x 10"' 3.98 1.248lQ;oJg l-365to ol8 0.085 0.3 

B1913+16 59.03 8.63 7.752 l-4408lo°°o3 l-3873loo°o3 0.617 0.3 

J1756-2251 28.46 1.02 7.67 l-40+o;°3 I-I8I003 0.181 1.7 

B1534+12 37.90 2.43 10.098 1.3332toool 1.3452lg;°°} 0.274 2.7 

J1811-1736 104.182 0.091 451.20 l-62lo;g5 l-Hto.is 0.828 >10 

J1518+4904 40.935 0.02 207.216 l-56l°;4g l-05^o.ii 0.249 >10 

J1829+2456 41.0098 0.05 28.0 l-l^lo.Ig l-36lo:i? 0.139 >10 



iGondek-Rosinska et al.l ()2005bl ) the authors suggest that 
there will be differences between the two observed popula- 
tions. In particular they reached the result that the mass ra- 
tio distributions will not be the same. The population mod- 
elling in those papers was quite crude. The binaries with 
lifetimes above the Hubble time were neglected, and it has 
been assumed that observability of pulsars is proportional 
to its lifetime. Such difference between the observed pop- 
ulation in the radio in gravitational waves would be quite 
important as we need to know how the gravitational waves 
signal looks like to detect it. In order to know this signal we 
need the initial parameters for numerical relativity calcula- 
tions. Until recently most calculations of coalescing neutron 
star binaries and black hole binaries assumed that mass ra - 
tio q is close to unity (e.g. iGondek-Rosihska et al.l (|20o3), 
iBuonanno et al.l l|2006h '). In this paper we present carefully 
model the binary p ulsar population with t he updated stel- 
lar evolution code (iBelczvnski et af] I2OO8I ) combined with 
the detailed pulsar evolution model. 



2 DESCRIPTION OF THE MODEL 

2.1 StarTrack binary population synthesis code 

In order to perform population synthesis of the binary pul- 
sars we needed to know how the population of NS-NS bina- 
ries itself looks like. We used th e StarTrack code for this 
purpose l|Belczvnski et af]|2008l ) . This code has been im- 
proved for many years now and is considered to be state- 
of-the-art with regard to our current understanding of stel- 
lar evolution. StarTrack has been extensively tested, cal- 
ibrat ed and its resul t s wer e compared with ob s ervati ons 
(e.g. 'Belczvnski et al' (I2OO2I ') . iBeTczvnski fc Buli^ (I2OO2I ') or 
[Bclczvnski ct al. (200^" ). At the base of this code are 
the analytic for mulae describing evo lution of single stars 
first developed bv lHurlev et all (|2000l ). The code follows se- 
quence of evolutionary phases of stars starting at main se- 
quence, passing through Hertzsprung gap, red giant branch 
or asymptotic giant branch. After the nuclear evolution of 
star ends it can form a compact object (white dwarf, neu- 
tron star or black hole). It can be also completely disrupted 
during supernova explosion leaving no remnant at all. Star- 
Track includes also the evolution of helium stars. If a star 



which first evolved ends up as a white dwarf or a neutron 
star then there is a chance for formation of a double neutron 
star system, provided that the binary is not disrupted during 
supernovae explosions. In the current version of StarTrack 
IBelczvnski et al]|2008l ) in a standard model we assume bi- 
modal distribution of the FeNi c ore masses in pre-supernova 
stars (Weaver & iWeaver fc Wooslevi (1993), T immes et al.l 
(| 19961 )') stars with initial masses below ~ 18 — 19 Mq burn 
carbon convectively forming cores of ~ 1.5 (NS of 1.3 M©), 
while stars of higher mass burn carbon radiatively and the 
cores are more massive ~ 2 (NS of 1.8 Mq). The mass of 
the FeNi core does not depend sensitively on the initial star 
mass, but is a strong function of burning reactions within the 
pre-supernova star, and if similar reactions are encountered 
over a wide initial star mass range in which neutron stars are 
formed (e.g. 8-18 Mq) then the majority of neutron stars are 
expected to have similar masses. The previously employed 
neutron star mass calculation was based on stellar evolution- 
ary modeling that was suggestive of in crease of final FeNi 
core m ass with initial star mass (e.g.. Weaver fc Wooslevi 
j 19931 ): or see Table 4.4 of IVanbeveren et al.l (|l998!)) ; that 
naturally resulted in a wide neutron star mass spectrum. 

Some poorly known details of stellar evolution are 
parametrized in the StarTrack code. Fortunately most of 
them do not have significant impact on the results relevant 
for this work. The various models of binary code used in this 
paper are described in section [2. 2. II 



2.2 Evolution of the binaries 

We consider three variants of the binary population synthe- 
sis models. They are described in detail below. 



2.2.1 Model A - standard 



This model is described in detail in IBelczvnski et al.l (|2002l) 
and updated in IBelczvnski et al.1 (|2008l ). Here we present 
only its most important features. Initial mass of heavier 
(at the ti me of form ation) star is drawn from power law 
~ M-^-^ IScaldl 19861 ). Mass of the second star is a fraction 
of it randomly chosen from a fiat distribution (|Kuiperlll935l ) , 
the initial orbital separation is flat in logarithm, and we draw 
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the initial eccentricity from a flat distribution. The code cal- 
culates the carbon-oxygen core mass which will form a com- 
pact object u sing formulae from already mentioned work by 
iHurlev et alj ( )2000. ) . This mass has impact on the mass of 
newly formed neutron star d ue to its influen c e on a mass 
of Fe-Ni core as described in JWooslevI (|l986h . iFrver et alJ 
l| 19991 ) and lTimmes et all (|l996l ) . 



MpeNi , 

M, 

Mf^n^ + fb{M ^ 



M 



FeNi), 



Mco ^ 5M0 
Mco > 7.6M0 
otherwise 



where ft is a fall-back factor (describing what fraction of 
the stellar envelope falls back), MpeNi is Fe-Ni core mass, 
Mco is CO core mass and M is pre-SN mass of the star. 
Orbital parameters are affected by tidal interactions, mass 
loss due to stellar winds, mass transfers, and gravitational 
waves emission. We parametrize the fraction of mass (/a = 
0.5) taking part in stable mass transfer that is accreted on 
companion and the rest is lost from the binary. Unstable 
mass transfer in StarTrack (common envelope - CE ) are 
treated with th e formalism proposed bv lWebbinkI (|l984l ) and 
Ide Kooll (I1990I ). In this model we assume acE x X = 1.0. 
Neutron star can accrete some mass during the common 
envelope phase. The amount of this mass is drawn from a 
range [O.O5M0,O.1OM0] with flat distribution. Supernova 
explosions are treated with details: the explosion occurs at 
a random pl ace in the orb i t, we take into account natal kicks 
drawn from iHobbs et al.l (I2005D distribution, and calculate 
the orbit and center of mass velo city after the explosion . 
This model is described in detail in iBelczvnski et "all (|2008l ) 
and is referred further as the model A, the standard model. 
The following two section describe two other models that 
differ from A. 



2.2.2 Model H 

The common envelope phase is treated differently in this 
case so that a neutron stars accrete more matter. It has 
crucial impact both on binary evolution (changes in or- 
bital separation) and on properties of the pulsar which 
formed first. In this m odel we assume Bondi-Hoyle accretion 
Bondi fc Hovld Il944l ) with a hypercritical accretion rate 



Brown et al 



]S) 



. The amount of accreted mass is roughly 
10 bigger than in the standard model A. This amount de- 
pends mostly on mass ratio in the binary prior to the com- 
mon envelope qp. Fin al mass is calculated using formula 
l|Bethe fc Brownll 19981 ): 

Mb = {qo + 0.7(?o) Ma 

where Mb is the mass of the acceptor after the transfer and 
Ma is the mass of the donor before the transfer. 



2.2.3 Model S 

In this model we assume a different initial mass function 
of newly formed neutron stars. It assumes very simplified 
linear relation between Mns and the mass of the iron core 
Mcore from which it is formed: 

Mns = 0.35 x M^ore + O.596M0. 

The values of parameters are such that we have neutron stars 
with masses from l.lAf0 (forming from a star with mass 



8.275M0 on a zero age main sequence) to 2.5M0 (form- 
ing from a star with mass 2O.88M0 on a zero age main se- 
quence) . The purpose of this model is to check the influence 
of this poorly known relation on double neutron star mass 
distributions. 



2.3 Pulsar evolution model 

In this section we describe the modelling of processes impor- 
tant in evolution of pulsars in binaries. At the very beginning 
we assign each neutron star a set of initial parameters such 
as spin period, magnetic field and moment of inertia. After 
that we start evolving each pulsar. First of all we need to 
know how its rotation behaves in time. We assume dipole 
model for the magnetic field. We also present how mass 
transfers influence neutron stars (e.g. recycling). During ac- 
cretion neutron stars are visible in X-rays but in this paper 
we neglect that as this periods are short compared to radio 
activity periods. We finish the evolution of a binary pulsar 
either at the time of coalescence due to gravitational wave 
emission or when it reaches the age of 10 Gyr. During the 
evolution we note the radio loud periods in the history of the 
binary, and also the periods it is visible in the radio from 
Earth. 



2. 3. 1 Initial parameters 

Each neutron star is evolved as a pulsar right after it 
is formed. Its initial magnetic field Bo is drawn from 
range [lO^^G, 10^'^G] with a fl at logarithmic distribution 
(|Faucher-Giguere k, Kaspill2006l ) . Each pulsar is ass igned an 
initial spin period Pini = 10 ms l|Cordes et al.|[20o3 ). We as- 
sume that each pulsar has the same Pint because even if 
some pulsars are born with longer or shorter periods then it 
does not change the overall properties of the population. It 
is so due to the fact that rotation of a young pulsar evolves 
very quickly so that a young neutron star on a short time 
scale moves on the P — P diagram towards higher spin pe- 
riods.. All the neutron stars are of equal radius J? = 10 km 
and a moment of inertia I = lO**^ g x cm^. We consider 
the spin-down only to to the electromagnetic radiation so 
the braking index n=3. The initial rotation frequency time 
derivative is: 



B^W sin^ a 
6c37 



(1) 



where a is the angle between rotation and magnetic field 
axis, B is the magnetic field, R is the neutron star radius, 
c is the speed of light and / is the moment of inertia of the 
neutron star. We assume a = 90° . for each pulsar. The mag- 
netic field value and the angle a are degenerate so we are 
not losing anything in making this assumption. The obser- 
vations show that t he breaking indexes of p ulsars are in the 
range ~ 2.5 — 3.5 (jManchester et al]|2005l ). At every point 
of evolution the magnetic field is given by 



B = 3.2 X 



2.3.2 Evolution of the magnetic field 



(2) 



There are no direct interactions between companions after 
formation of both neutron stars (except for coalescence and 



Population synthesis of binary pulsars. 



5 



contact of magneto-spheres). At this stage the spin period 
and its evolution do not depend in any way on companion's 
properties. If this statement is true for a given numerical step 
then we call it a standard step. We can have such steps also 
in the time between formation of first and second neutron 
star if there are no mass transfers ongoing. In the standard 
step we calculate current spin period and its derivative ac- 
cording to equation [T] In addition we assume magnetic field 
of a pulsar is decaying exponentially in time: 



radius equals: 



B = {Bo — B„ 



i) • exp + B„^i„ , 



(3) 



where Td is the magnetic field decay time scale. The 
decay of the magnetic field is controversial with some 
analyses sho wing that it p r oceed s on a timescale of 
2 - 5Myrs (iGonthier et all |2002| . 120041 ') . while other 



iFaucher-Giguere fc Kaspi 20061 ) arguing that it may be an 
artifact of the assumed luminosity law. We include it as in 
general there seems to be some torque decay in pulsar evo- 
lution (jFaucher-Giguere fc Kaspill2006l ). The pulsar which 
formed earlier can be an exception as we assume that after 
it accretes any amount of mass its field ceases to decay. We 
also assume that there is a lower limit for the strength of 
magnetic field Bmin = 10* G., i.e. that the magnetic field 
decay ceases once the field reaches this value. 

2.3.3 Roche lobe overflow 

Before the second supernova explosion occurs, the compan- 
ion of an already formed pulsar continues its stellar nuclear 
evolution. It can fill its Roche lobe in this period. If that 
is the case then there is mass transfer from the companion 
and the system becomes an X-ray binary. Some mass is lost 
from the system and some is accreted onto the neutron star. 
Besides increase of the mass of the neutron star mass, the 
mass transfer results in three effects: 

(I) we assume that the magnetic field decays, and we model 
the magnetic field change as: 

B = {Bo - Bmin) ■ exp ^--^^^ + Bmin , (4) 

where dM is the accreted mass and AMd is magnetic field 
decay mass scale , and Bmin is the minimal magnetic field of 
a neutron star. After the mass accretion induced decay the 
field stops to decrease in time, i.e. it is fixed at the final level 
unless another mass transfer happens. In the latter case field 
decays in the same way as in the first event. The lower limit 
for the magnetic field also applies Bmin ~ 10* G. 

(II) the pulsar in the accretion phase is not emitting in the 
radio. This is because the magnetosphere is filled by the 
accreted gas, and the particle acceleration is not effective. 
The pulsar is visible in the X-ray band but in this work we 
do not model this phase. 

(III) the spin period of a pulsar is affected. Due to angular 
momentum transfer the neutron star tends to corotate with 
the accreted matter at the Alfven radius: 



Ra 



/8R^y 
\MAPG J 



The accretion rate M is calculated within the StarTrack 
code by considering the de tailed evolution of the compan- 
ion l|Belczvnski et alll2008l ). The orbital velocity at Alfven 



0/ 



GM 



If the amount of matter is sufficient, namely dM > Q.IMq, 
then we assume that such corotation is attained and the 
pulsar is fully spun up. Note that if the matter at Alfven 
radius has angular velocity smaller than this of the neutron 
star, it can spin down the neutron star. In such cases the 
propeller effect might take place (see 12.3.6"]) . 

2.3.4 Partial spin-up 

If the accreted mass is smaller than O.IM© then the pulsar 
can be spun-up to a smaller angular velocity. We assume 
that the amount angular is momentum transferred AJ is 
proportional to the accreted mass: 

AM 
A J cx , 

O.lMg ' 

This leads to a linear relation between the amount of ac- 
creted mass and the final angular velocity of the pulsar: 

AM 



Q.f = 0.0 + {Q.A - fio); 



O.IMq ' 

where f2/ is the final angular velocity, Q,a is the angular 
velocity at the Alfven radius and Oq is the initial angular 
velocity before recycling. 

2.3.5 Common envelope 

Formation of double neutron stars involves quite frequently 
an unstable mass transfer, i.e. a common envelope phase. 
During the common envelope phase the binary orbit is tight- 
ened, and a significant amount of matter from the com- 
panion is expelled. At the same time some matter is ac- 
creted onto the neutron star, see section 12.2.11 and refer- 
ences therein. In the St arTrack code th e common envelope 
is treated following the (|Webbinklll984l ) formalism and it is 
assumed that the amount of matter accreted by the neutron 
star lies in range [O.OSAfg , CIOA^q] and is drawn from it 
with a fiat distribution, which corresponds to highly super 
Eddington accretion rat e. At this stage the ang ular momen- 
tum transfer is chaotic (|Benensohn et allllQQTl ). and we as- 
sume that the spin of a neutron star is not affected on the 
average by accretion during this phase, i.e. the spin period 
of the pulsar is the same as before the common envelope. 
However, the accreted matter quenches the magnetic field 
of the neutron star, according to equation (Q. Thus, a sys- 
tem after a common envelope phase has a tightened orbit, 
and the pulsar is more massive with decreased magnetic field 
yet its spin frequency is unchanged. 

2. 3. 6 The propeller effect 

Here we consider the case of a pulsar rotating with 
greater velocity than matter orbiting at the Alfven radius. 
llUarionov fc SunvaevI l|l975h suggested that in this case the 
propeller effect takes place: accretion onto neutron star is 
inhibited because of centrifugal force at the Alfven ra dius. 
This hypothesis is supported by the X-ray observations (|Cuil 
Il997f) . 
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In section [2.3.31 we describe how the accretion can spin 
up or down the neutron star. In models which include the 
propeller effect, we do not allow for the spin-down by ac- 
cretion, however, we allow for the magnetic field decay and 
with no infiuence on the spin rate of the neutron star. 

In the models where we neglect the propeller effect a 
neutron star can both be spun down or up. The outcome 
depends on the value of the magnetic field and on the ac- 
cretion rate. 



2. 3. 7 Evolution of the orbit 

Some of the binaries are tight enough to be strongly in- 
fluenced by gravitational waves emission. We calculate the 
orbit evolution due to emission of gravitational waves in first 
Post-Newtonian approximation. The time derivatives of the 
eccentricity and the semi-major axis are (Pet ers 19641 '): 



da 
IE 



64 G'\Nh + M2f 
5 c'a^il — e^)i 



l+e' 



/73 37 2\ 
(24 + 96^) 



de 
IE 



304 G-HM-i + Mi 
-e 



f, ^ 121 2 
1^+ 304^ 



5 c5a4(l-e2. 

Orbit evolution is taken into account for the entire time the 
binary is evolved, namely for 10 Gyr or until coalescence. 
The pulsars can be radio loud until the moment when the 
semi-major axis of system is smaller than the sum of light 
cylinder radii: 



a < 



where a is the semi-major axis, c is the speed of light, Qi 
and Q2 are angular velocities of first and second pulsar re- 
spectively. The contact of magneto-spheres in practice does 
always happen in the same numerical step as the coales- 
cence. 



2.3.8 Radio luminosity 

In order to check the infiuence of selection effects and to 
compare populations visible in gravitational waves and in 
radio we n eed to know the radio luminosity of all the pulsars. 
We use a iNaravan fc Ostrikeij l|l990D model fitted to the 
observations: 



logio L = - logio 



p3 



+ 1.635, 



where P-15 



This luminosity is given in units of 



mJy X kpc^ for observations at 400 MHz. 

The pulsars cease to emit in the radio once they cross 
the death lines. After crossing this lines on this diagram the 
pulsar's emission mechanism fails as electron-positron pairs 
can no longer be create d in the magnetic field . We take two 
death lines into account (|Rudak fc Ritterll994h . In this work 
we assume that all the pulsars cease their emission after 
crossing the death lines defined by: 

logio P = 3.29 X logio P - 16.55 , 

logio P = 0.92 X logio P - 18.65 . 
These relations are rather empirical and there are cases when 



a pu lsar is found beyond these death lines (|Young et al.l 
Il999h . yet they describe the cutoffs on the P-P diagrams 
quite well. 



2.3.9 Numerical step size 

In order to choose the length of the numerical step size dur- 
ing the evolution of a DNS we calculate a number of time 
derivatives describing the rate of change of pulsar proper- 
ties. The step is chosen in such a way to satisfy the following 
conditions: (i) the pulsar period changes by less than 1% (ii) 
the semi-major axis of the orbit changes by less than a 1%, 
(iii) the step is limited by the beginning of the next mass 
transfer episode, (iv) the time step is shorter than 1 Gyr. 
The last condition was introduced to follow the evolution of 
recycled millisecond pulsars with weak field on large orbits. 



2.3.10 Free parameters 

We list the parameters describing the evolution of pulsars 
in Table |2l We consider several models of pulsar evolution 
and we list them along with the values of the evolutionary 
parameters in Table O We choose the following notation to 
denote each model: the first letter in the name of model 
refers to the stellar evolution model in the StarTrack code; 
the letter P means that the model incorporates the propeller 
effect. The letter F implies that no partial spin-up is possible 
in given model, all the pulsars are recycled to the orbital ve- 
locity at Alfven radius regardless of the amount of accreted 
matter; the letters D and T are followed by numbers refer- 
ring to the magnetic field decay - respectively mass(AAfd) 
and time (rd) scales. 



2.3.11 Motion in the Galaxy 

In order to obtain the observed radio fiuxes of the pulsars 
we need to model their motion in the Milky Way. Therefore 
we consider a simple model of the Galaxy consisting of the 
following three components: bulge, disk, and halo. The bulge 
and d isk potential are described by th e Mivamoto fc Nagail 
(|l975l l type potential (jPaczvnskilll990l : iBulik et al.lll999l ): 



$(r, z) 



GM 



y/R^ + {a + Vz'' +b''y 



where M is the mass of a given component, R = y^x^+y^, 
and a, b are the parameters. The halo is described by the 
density distribution p — pc[l + {r /rc) with a cutoff at 
rent = 100 kpc above which the halo density is zero. The 
corresponding potential for r < rcut is 



$(r) 



GMh 



1 , / r^\ rc 

— In 1 H — 7; H arctan 

2 V ri r 



We use the tol l owmg values ot tiie parameters 
l|Blaes fc Raiagopall Il99ll ) describing the bulge (index 
1) and disk (index 2) potential: a\ — Okpc, bi = 0.277 kpc, 
a2 = 4.2kpc, b2 = 0.198kpc, Mi = 1.12 x 10^" Mq, 
M2 = 8.78 X 10^" Mq, while for the halo potential we use 
Mh = 5.0 X 10"' Mq, and = 6.0 kpc. The distribution 
of stars in the disk is assumed to be that of a young disc 
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Table 2. Free parameters in the evolution model 



Parameter Description Default value 

r^j Time scale of field decay 5 Myr 

AMd Mass scale of field decay 0.025 Mq 

Bo Initial magnetic field [lO" G, 10^^ G] 

Bmin Minimal magnetic field 10* G 

/ Moment of inertia 10*^ g X cm'^ 

R Pulsar radius 10 km 

a Angle between rotation and magnetic axis 90° 

— Mass limit for full spin-up 0.1 M© 



l|Paczvnskilll990l ). The radial and vertical distributions are 
independent i.e. the distributions factor out: 

P{R, z) oc R{R)dRp{z)dz , 

where the radial distribution is p{R) oc exp{—R/ Rexp), 
and Rexp ~ 4.5 kpc, and we introduce an upper cutoff at 
Rmax = 20 kpc. The vertical distribution is exponential 
p{z) oc exp (— z/75pc). After distributing the binaries with 
the given densities, we evolve them in time. We use the 
leapfrog method in KDK scheme with a constant (10^ years) 
time step. This approach is fully symplectic. 

2.4 Detection in the radio band 

We consider radio population to be composed of pulsars 
with non-zero radio luminosity. The radio luminosity of each 
pulsar is calculated at each time step. The radio luminos- 
ity of a pulsar drops to zero during mass transfers and af- 
ter crossing the death lines. We define the radio population 
with selection effects as all objects with the ffuxes at Earth 
above ImJy. Additionally the BNS with the orbital pe- 
riod between .3fe and 4/i may have decreased detectability 
l|Faulknerl2004l ). We consider a population with and without 
such objects and we will discuss the differences (see ll.2|l . We 
neglect the fact that radio beams can have a limited width 
and the dispersion by electron gas in the Galaxy. 

2.5 Detection in gravitational waves 

The gravitational wave detectors su ch as LIGO: 
(.Fla nagan fc Hughes, ,1998) and VIRGO (|Kulczvcki et al.1 
l2006l ) are currently working and in principle can detect 
neutron star coalescence as far as 18 Mpc. The signal to 
noise ratio from the inspiral phase of two compact objects 
is given by 

— ~ X — 

where D is a distance from the observer and chirp mass 
Mchirp = (MiM2)°'® (Ml -I- M2)"°'^, where Mi and M2 are 
masses of first and second star in the system respectively. 
We define the DNS population observable in gravitational 
waves as these with the merger time below the Hubble time 
and weight all their observable quantities with the volume 
in which they are observable, i.e. 

v~m|- . 

chirp 

Thus the gravitational wave population is the population 
residing in multiple galaxies and we assume that the pop- 
ulations in these galaxies resembles the one in the Galaxy. 




-3 -2 -10 1 



log(P) 

Figure 2. An example of a neutron star binary evolution on the 
P ~ P diagram. The dotted line corresponds to the birth line, 
where all neutron star are formed in our model. The solid lines 
represent the death lines. The numbers correspond to character- 
istic ages in Myrs and are explained in in the text 1)3.1^ . Arrows 
show how the mass transfer affects the first pulsar. 

We neglect the possible detection of the gravitational waves 
in the merger and ringdown phases of the coalescence. 



3 RESULTS 

3.1 Example of a binary evolution 

In subsection 12. 31 we described our phenomenological model 
of pulsar evolution. In Fig. [2] we present an example of an 
evolutionary path leading to the formation of a double neu- 
tron star. The example is based on the model AF. 

The pulsars are born with the rotational period Pini — 
10 ms at the line of birth (the dotted line on Fig. [2}. The ini- 
tial value of the magnetic field ,which determines the spin 
period derivative, is drawn from a fiat distribution as de- 
scribed in section 12.31 The binary starts on zero age main 
sequence at t = 0. The first pulsar is born after 28.08 Myrs. 
The system consists of a pulsar and a massive rejuvenated 
companion. The neutron star evolves in the P — P plane 
along the line of constant magnetic field. At t=31.26 Myrs 
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Table 3. Models of pulsar evolution 



Model 




AMa 


propeller 


spin-up 


StarTrack model 


A 


5 


0.025 


no 


partial possible 


A 


AF 


5 


0.025 


no 


only full 


A 


APD05 


5 


0.05 


yes 


partial possible 


A 


AP 


5 


0.025 


yes 


partial possible 


A 


APF 


5 


0.025 


yes 


only full 


A 


APT20 


20 


0.025 


yes 


partial possible 


A 


HP 


5 


0.025 


yes 


partial possible 


H 


SP 


5 


0.025 


yes 


partial possible 


S 



the nuclear evolution of the companion plunges the system 
into the common envelope phase. The magnetic field of the 
pulsar is quenched and it falls close to the death line. The 
companion loses the envelope and becomes a helium star 
while the neutron star reappears in radio. It barely evolves 
until t = 37.11 Myrs when the second mass transfer occurs, 
this time it is a stable Roche lobe overflow. The system be- 
comes an X-ray binary and the neutron star is recycled to a 
period ~ 0.01s. The amount of accreted matter is very small 
(~ 10~^Mq) but in this model we allow for full recycling 
anyway. In other models we limit the accretion rate to the 
Eddington limit and consequently the recycling is weaker. 
Soon after this mass transfer the second pulsar is born, and 
for the following 23 Myrs the binary contains two radio loud 
pulsars. At around t=50.74 Myrs the second born pulsar 
falls below the death line. At this stage the system contains 
the recycled pulsar which slowly evolves to the death line. 
It does not reach it throughout the whole simulation, which 
ends after 10 Gyrs have passed. Some of the pulsars in our 
simulations manage to pass the death line before coales- 
cence. This binary is relatively wide, so it does not coalesce 
within that time. 

3.2 Standard model 

In order to analyse the properties of the population of binary 
pulsars we have to take into account the selection effects. We 
define two types of selection effects that leads to differences 
between the intrinsic and the observed population. The first 
is the simple effects connected with the fact that any radio 
survey has limiting flux. We select the radio population as 
all pulsars with the observed flux on Earth F > ImJy. Ad- 
ditionally, it has been noted that binaries with the orbital 
period between 0.3 and 4h are difficult to detect (Faulkner, 
2004) so we also remove these systems from the observed 
population unless otherwise noted. 

3.2.1 Population on P — P 

We start with the analysis of the population of binary pul- 
sars on the P— P diagram. To obtain the present observable 
population of binary neutron stars we assume that star for- 
mation rate in the galaxy is constant. We place them in the 
Galaxy according to the model described in section 12.3.111 
and we propagate them in the Galactic potential taking 
into account the additional kick velocities that the binary 
receives at birth of the two pulsars. We then look for the 
observable population taking into account the selection ef- 
fects described above. The density of the objects in the P~P 



plane is found by summing the times that each pulsar spends 
in a given cell on the log P — log P plane, i.e. the density in 
the in a given cell on P — P is : 

F 1 

/ ^ £21 i , , (5) 

■' ' F AlogPAlogP ^ ^ 

where AlogP = 0.08 and AlogP = 0.18 denote bin width 

50 

of spin period and its derivative respectively, P — X/ 
and 

where Nij is a count of simulated pulsars with P contained 
in the i-th and P in j-th bin and tf^ is a time p-th pulsar 
spends in these bins. Range of logj^Q P from —3 to 1 as well 
as logj^Q P range from —21 to —12 were divided into 50 bins. 

The density of the population in the P—P plane (shown 
on Fig. [3] and |3| is determined by the properties of the two 
component of the population: the first born, possibly recy- 
cled pulsar and in a much smaller degree by the second born 
pulsar. The details of the binary evolution and recycling are 
influencing strongly the region where we expect the recycled 
pulsars. For the models where we allow for the full spin-up 
(AF, APF) the population of millisecond pulsars extends to 
period below 20ms. Another model where there are many 
pulsars in this region is model HP, where we allow for very 
strong accretion and therefore recycling is always full. In 
model HP most pulsars end up in the small region with 
the pulsar period of approx l-2ms, and period derivative of 
~ 10~^^. In the remaining models the mass transfer leads 
to decrease of the field in the CE phase, but at the same 
time the period is not changed by much. Inclusion of the 
propeller effect (models denoted by P) limits the amount of 
accreted matter and thus may limits the field decay. This 
decreases slightly the density of pulsars with the lowest val- 
ues of P, see e.g. model A vs model AP. The time scale of 
the spontaneous decay of magnetic field does not influence 
the overall shape of the distribution in the P — P diagram 
see model AP and APT20. Finally, the magnetic field mass 
decay scale influences the values of the magnetic flelds of the 
recycled pulsars and hence the distribution of their period 
derivatives P, see model APD05 with the AM^ = 0.05 Mq 
and AP with AM^ = 0.025 M©. If the value of AMd is 
large than the magnetic field of pulsar does not decay fast 
and they end up as millisecond recycled pulsars with a lit- 
tle stronger magnetic field in the region between 10~^^ and 
10~^^, while for the small value of AA/d the final magnetic 
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Figure 3. Probability density on P — P diagram within models APD05,AF,APF, and APT20. The levels correspond to contours 
containing 68, 95, and 99 percent of the objects. The outer contour delimits the region containing all the objects in the simulation. Black 
dots correspond to observations. 



field of pulsars are lower and they populate the region below 



10" 



see model AP. 



3.3 Comparison with observations 

In order to quantify which model of the distribution in the 
P — P best describes the data we calculate the likelihood 
of each model given the data from Table 1. We define the 
likelihood as : 



where C is the likelihood, Pj, and Pk are period and its 
derivative values for a k-th observed pulsar (see Tab.[T} and 
/ is the probability density defined by equation [5] in section 
13.2.11 We use two sets of data to calculate the likelihood: 
(Comparison 1) the standard model, where we take into ac- 
count the flux selection effect as well as rejection of pulsars 
with the orbital periods between 0.3 and 4h is compared with 
the data set with J0737-3039 A and B excluded as its or- 
bital period lies in this range; (Comparison 2) the extended 
model where we do not reject the pulsars with the orbital 
periods between 0.3 and 4h is compared to the full data set 
of Table [T] and (Comparison 3) with the data set where we 
exclude J0737-3039B. We present the results in Table g] 
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-1 

log(P) 

Figure 4. Same as Figure [S] for the models AP,A, SP, and HP. The yellow contour for HP model overlaps with the green one. 



In all three comparisons in Table [53] the model APD05 
best fits the data. This is a model with the propeller effect 
taken into account and with increased mass scale for the 
magnetic field decay. The models with the standard value 
lead to too small magnetic field of the synthetic sample. 
Increasing the magnetic field mass decay scale AMd leads 
to smaller decay of the field, and therefore better agreement 
with the data. Also the models with full recycling are very 
far from the observations since they predict the population 
of pulsars in binaries with spin periods below 20ms. 



EXPECTED MASSES OF NEUTRON STARS 
OBSERVED IN GRAVITATIONAL WAVES 
AND IN THE RADIO 



In iBulik et al] l|2004l ) the authors have calculated the ex- 
pected masses of neutron stars observed in gravitational 
waves using the StarTrack population synthesis code (Bel- 
czynski et al. 2002 for detailed description). The assumed 
value for the minimum mass and maximum mass of a neu- 
tron star was 1.2 Mq and 3.0 Mq respectively. They have 
verified the res ults taking the maximum mas s to be 2 Mq 
and 2.5 Mq l|Gondek-Rosiiiska etahl l2005ah . They have 
shown that the distributions of mass ratio q (defined as the 
ratio of less to more massive component) of neutron star 
binaries observed in gravitational waves, have two peaks: 




Figure 5. Distribution of masses (in Mq) and mass ratios observed for the radio (left panels) and gravitational wave samples (right 
panels). The shadings represents the region containing 68% , 95% , and all the systems. In each panel we show the lines corresponding 
to constant chirp masses of 1.0, 1.1, 1.2 A/q. Black dots correspond to observations. 
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Table 4. The logarithm of the hkehhood for all models. The 
missing entries correspond to to the case where the likelihood is 
formally — oo due to the fact that we find no model pulsars in the 
bin containing J0737-3039B. 



Comparison 


1 


2 


3 


A 


-14.24 


-20.31 


-16.97 


AF 


-19.90 




-23.31 


AP 


-16.52 




-18.27 


APD05 


-14.21 


-19.26 


-16.00 


APF 


-21.13 




-23.06 


APT20 


-15.52 


-21.27 


-17.65 


HP 


-22.28 


-27.36 


-24.14 


SP 


-15.27 


-20.11 


-17.26 



the first one for nearly equal mass systems with two masses 
close to 1.4 Mq; and the second one with small mass ra- 
tio q ~ 0.6 — 0.7 with binaries consisting of a neutron star 
with gravitational mass close to the maximum mass with 
relatively smaller companion. The strength of the second 
peak depends on the value of the maximum mass of a neu- 
tron star. For the models with the higher maximal mass the 
sample starts to be dominated by binaries belonging to the 
second peak. The authors used 20 extra different models of 
stellar evolution, where they have varied the parameters de- 
scribing various stages of stellar and binary evolutions in 
order to assess the robustness of the results. The binaries 
with a small mass ratio has so far been overlooked in rel- 
ativistic calculations of the hydrodynamical inspiral or the 
merger phase. In that paper, the gravitational wave selected 
sample was compared with the volume limited sample. The 
volume limited sample contained only the binaries poten- 
tially detectable in gravitational waves, i.e. those younger 
that lOGyr. 

Up to now there was no detailed comparison of proper- 
ties of neutron star binaries observ ed in radio and in grav- 
itation al waves. The only paper l|Gondek-Rosiiiska et al.l 
l2005bl ) that dealt with a similar problem contained a discus- 
sion of radio observability of the gravitational wave selected 
sample, within a single model of binary evolution with a wide 
initial neutron star mass spectrum. The authors assumed 
that the radio selected sample consists only of binaries with 
lifetimes (from formation till merger) longer than 100 Myrs, 
and there was no upper limit on the radio lifetime other 
that the Hubble time. They were observable as pulsars for 
the entire lifetime of the binary. In that paper it was shown 
that binary radio pulsars, which are extended systems, were 
only a few percent of all binary neutron stars observed in 
gravitational waves. The binary pulsars observed in the ra- 
dio had mass ratio distribution clustered around unity. The 
long lived systems evolved without the possibility of signif- 
icant accretion onto a neutron star while the short lived 
systems (with tgrav < 100 Myrs) did undergo common enve- 
lope episodes with hypercritical accretion onto the neutron 
star. This common envelope episodes had two consequences: 
they tightened the orbits and led to decrease of the mass ra- 
tio of the final double neutron star system, as one of the 
neutron stars accreted some matter. This implied that the 
mass distribution of the gravitational wave selected sample 
of double neutron star binaries was different than the radio 
selected one. 



Here we present a comparison of the radio selected sam- 
ple of double neutron star binaries, taking into account selec- 
tion effects, with the gravitational wave selected one. The ra- 
dio population contains also the binaries with merger times 
above the Hubble time. 

In FigureOwe present the distributions of expected ob- 
jects in the plane spanned by the mass ratio and the primary 
mass defined as the more massive component of a binary for 
APD05, HP and SP models. The left panels correspond to 
the radio selected sample while the right panels show the 
gravitational wave selected ones. The objects contained in 
the radio sample are weighted by the time they are observ- 
able as pulsar from Earth, while the objects in the gravita- 
tional wave selected sample are weighted by the volume in 
which they are observable. The solid lines correspond to con- 
stant values of the chirp mass in these coordinates. Observed 
binary neutron stars are shown as black dots. In addition on 
Figure |S] we show the chirp mass distributions of the bina- 
ries selected by their observability in the radio band, and by 
their observability in gravitational waves. 

Comparing our radio selected sample with the data from 
Table 1 we see that the APD05 (top panel) and SP (lower 
panel) models are more consistent with observations than 
the HP (middle panel) model in which majority of binary 
radio pulsars are predicted to have g < 0.8 and one massive 
neutron star with m2 > 1.6Mq. The APD05 model, which 
reproduces best the observed distribution of pulsars in the 
P — P diagram, is consistent with most of the observed bi- 
nary radio pulsars. It does also reproduce very well the chirp 
mass distribution of the radio observed sample. The binary 
neutron stars with low mass ratios and moderate masses are 
not predicted in this model. However, one must take into ac- 
count the fact that measurements of masses of J1811-1736 
and J1518-I-4904 carry quite large error bars. The SP model 
reproduces the distribution of masses and mass ratios of 
all observed neutron star binaries. In this model the dis- 
tribution of masses is wide. It predicts existence of binaries 
containing massive neutron stars, however the radio selected 
distribution is concentrated around q ~ 0.8 — 0.9 and pri- 
mary masses ~ I.SSMq. 

A comparison of radio selected sample with the gravi- 
tational wave one shows visible differences for SP and HP 
models and negligible for the APD05 model (see also Figure 
[6]). In the model APD05, the two samples are very simi- 
lar as expected. This is due to the fact that in the binary 
evolution model A the range of masses of newborn neutron 
stars is narrow comparing to the previous calculations as 
well the amount of matter accreted during the common en- 
velope episode with a helium star is negligible. 

In model HP the radio sample is dominated by unequal 
mass binaries with q — 0.7, m2 ~ 1.7Af0. However the grav- 
itational wave selected sample the equal mass binaries with 
both components of about 1.4M0 are also dominant. Also 
the chirp mass distributions are very different, with radio 
sample leaning towards the higher chirp mass binaries. This 
result which is apparently counterintuitive, follows from the 
fact that in this model we allow for substantial accretion and 
consequently strong recycling of the first born pulsar. This 
pulsar has low magnetic field and remains radio loud for a 
very long time, and thus it contributes significantly to the 
radio sample. At the same time the binary has a large chirp 
mass because of the strong accretion. The gravitational wave 
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selected sample is dominated in model HP by the short lived 
systems with no significant accretion episodes, and therefore 
contains more systems with a low chirp mass. Additionally 
this model fails to reproduce the radio observed distribution 
of masses and mass ratios. 

In the SP model the gravitational wave selected sample 
is shifted toward the higher chirp mass values with respect to 
the radio selected one. This is because of the volume effect, 
due to the fact that the sampling volume scales as oc M^^^. 
Thus the heavier binaries are observable from much larger 
volume in gravitational waves. The mass ratio distribution 
of the gravitational wave selected sample leans toward lower 
values because the unequal mass binaries typically contain 
a low mass neutron star with a more massive companion. 
Thus they have a higher chirp mass than the equal mass 
neutron star binaries that typically contain two stars with 
low masses. This is similar to the comparison between the 
gravitational wave selected and the volume limited sample 
presented in Bulik et al. 2004. The comparison with radio 
observations reveals quite and interesting feature: the distri- 
bution of chirp mass is underestimated in this model, how- 
ever the observed mass ratios are relatively well reproduced, 
as this model leads to a wide distribution of radio observed 
mass ratios. 

The small differences between the gravitational wave 
selected and the radio selected distributions of mass ratio 
and chirp mass in model APD05 are mainly due to the fact 
that the initial masses of neutron stars come from a very 
narrow range and we do not allow for substantial accretion. 
Thus most neutron star binaries have similar masses and 
there is no possibility to form binaries with different chirp 
masses or mass ratios in this model. However, once we allow 
for a wide distribution of neutron star initial masses or for a 
significant accretion the differences between the two samples 
start to appear. 



5 CONCLUSIONS 

We have modelled the evolution of binary stars leading to 
formation of binary neutron star. We then follow their evo- 
lution as pulsars and discuss the properties of the binary 
pulsars observable in radio. We discuss several models and 
compare the expected distribution with observations in the 
P-P diagram. 

First of all the analysis of the binary star evolution pro- 
vides little constraints on the spontaneous magnetic field 
decay. This was to be expected as binary pulsars are typ- 
ically observed already as old and evolved objects. At the 
same time we see that there must be some mechanism at 
work which leads to decay of the magnetic field of the first 
born star. This can either be the spontaneous magnetic field 
decay, mentioned above, or the magnetic field decay due to 
accretion. However we noticed that in the best fitting mod- 
els the amount of accreted matter is small. This is due to 
two factors: in some objects the propeller effect inhibits ac- 
cretion, while in other the time for the accretion is small 
because of the fast evolutionary timescale of the compan- 
ion. Moreover, the model with large mass decay scale fits 
the data better than the one with small value. Thus in bi- 
nary pulsars the spontaneous field decay seems to be much 
more important than the accretion induced decay. 
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Figure 6. The distributions of chirp masses for gravitational 
wave selected sample and for the radio selected sample of double 
neutron star binaries for model APD05 (top panel), HP (middle 
panel) and SP (bottom panel). 



The models in which we allow for substantial accre- 
tion and recycling of the neutron star lead to formation of 
a large number of binaries with periods below 20ms and 
period derivatives below 10~^*ss~^. There are no observed 
binary pulsars in this region which means that such substan- 
tial accretion does not take place. Thus the mass accreted 
by the first born star must be small, or the efficiency of re- 
cycling in common envelope is very low. This efficiency can 
be much smaller than in the case of stable mass transfer 
because the mass transfer can be turbulent and we do not 
expect a steady transfer of angular momentum. 

If the angular momentum transfer is low than we do 
not expect binary pulsars to have pulsation periods below 
~ 20ms. The binaries containing pulsars with lower period 
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arc probably mainly neutron star - white dwarf binaries, in 
which the mass transfer had been stable and lasted much 
longer than in the progenitors of double neutron star bina- 
ries. 

The orbital parameters of observed objects shows that 
the expected number of binary pulsars with orbital periods 
below 2 hours is very small. Thus, it is quite likely that the 
binary pulsar J0737-3039 with Port = 2.454 h will remain 
the lowest orbital period binary pulsar for a long time. On 
the other hand the models hardly can explain the orbital 
properties of J1518+4904. Its orbital period is typical for 
DNS with the eccentricity of w 0.8, much higher than the 
observed value 0.249. This eccentricity is also too large for 
a pulsar with a massive white dwarf companion. 

Finally, we compare the properties of the radio selected 
and gravitational wave selected sample. We find that two 
different models APD05 and SP are consistent with obser- 
vations. The model APD05 that best fits the distribution on 
the P — P diagram also reproduces well the distribution of 
chirp masses. In this model the binary evolution proceeds 
within the standard StarTrack model A, we allow for the 
propeller effect, partial spin-up proportional to the amount 
of mass accreted, and mass scale of magnetic field decay is 
increased to Md = 0.05 M©. In this model there are very 
little differences between the distribution of masses in the 
radio and gravitational wave selected samples, and the mass 
ratio of binary pulsars is close to unity in both the radio 
and gravitational wave selected sample. However, radio ob- 
servations show that there are systems with the mass ratio 
likely to be in the range between 0.6 — 0.9, which do not 
appear in this model. The presence of such objects can be 
explained with the assumption that the initial mass spec- 
trum of neutron star is wide like the one in model SP. The 
presence of such binaries suggests that also the unequal mass 
systems with the mass ratio in the range from 0.7 to unity 
should be included in the search for gravitational wave sig- 
nals. Inclusion of a wider initial mass range on one hand 
widens the expected distribution of mass ratios but on the 
other hand decreases the agreement between the expected 
and observed distribution of chirp masses in the radio. This 
conclusions may be verified with the Square Kilometer Ar- 
ray observations that should reveal many more binary NS 
systems. The measurement of their masses along with iden- 
tification of binary NS merger population by the Advanced 
LIGO and VIRGO will provide further insight into the prop- 
erties of the population of binary neutron stars. 
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